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Angular momentum projected multi-cranked 
configuration mixing for reliable calculation of 
high-spin rotational bands 


Mitsuhiro Shimada, Shingo Tagami and Yoshifumi R. Shimizu 


By employing the angular momentum projection technique we propose a method to reliably 
calculate the quantum spectrum of nuclear collective rotation. The method utilizes several 
cranked mean-field states with different rotational frequencies and they are superposed in the 
sense of the configuration mixing or the generator coordinate method, after performing the 
projection; the idea was originally suggested by Peierls-Thouless in 1962. It is found that the 
spectrum as a result of the configuration mixing does not essentially depend on chosen sets 
of cranking frequencies if the number of mean-field states utilized in the mixing is larger than 
a certain small value. We apply this method to three examples employing the Gogny DIS 
effective interaction and show that it is useful to study high-spin rotational bands by means of 
the angular momentum projection method. 
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1 Introduction 


The rotational motion in deformed nucleus is a typical symmetry-restoring collective 
motion and it appears as the rotational band, i.e., a group of eigenstates generated from 
one deformed intrinsic state [l|. In fact a nice rotational spectrum with I'" = O'*", 2~^, ■ ■ ■, 

can be obtained for the ground state band of even-even nuclei by the method of angular 
momentum projection from an axially-symmetric deformed state. The angular momentum 
projection is a standard quantum mechanical technique to restore the rotational invari¬ 
ance broken in the mean-held approximation, and is suitable to microscopically describe the 
nuclear collective rotation . However, this technique has not been very commonly utilized 
especially for the description of the high-spin states except for the works by the Projected 
Shell Model approach, see e.g. Ref. j^. In the present work we propose a method to reliably 
calculate the high-spin rotational band based on the angular momentum projection method. 

A standard method to describe the high-spin states is the mean-held approximation com¬ 
bined with the cranking model namely the Hamiltonian transformed into the uniformly 
rotating frame (the so-called Routhian), 


H — H l^iot Jyi 


( 1 ) 


is considered instead of the original Hamiltonian H. Here the quantity Wrot is the rotational 
(or cranking) frequency and Jy is the angular momentum operator around the rotation-axis 
(here chosen to be the y-axis). It was hrst introduced by Inglis to calculate the moment 
of inertia of the ground state rotational band treating the cranking term —UrotJy as a 
lowest order perturbation, i.e., at the Wrot —t 0 limit. It was extended to the hnite frequen^ to 
understand various phenomena at high-spin states as the selfconsistent cranking model or 
the cranked shell model j^, se e e.g . Refs. Recently its three-dimensional version, i.e., the 
so-called tilted-axis-cranking [l0|, was also successfully applied, see e.g. Ref. [ll|. The basic 
idea is that the selfconsistent mean-held state obtained by minimizing the Routhian in Eq. (IT]) 
depends on the rotational frequency, |*hcr(i^rot)), and is interpreted as an intrinsic state 
suitably describing the high-spin state at the spin value Jh ~ (*hcr(a;rot)|<Ti/|‘hcr(<^rot))- In this 
way the effects of the collective rotation on the mean-held parameters like the deformations 
and the pairing gaps can be well taken into account. 

The angular momentum projection is a fully microscopic method to calculate a sequence 
of quantum states belonging to one rotational band. Although we can obtain a whole sequence 
of the band by the projection from one intrinsic state, the resultant spectrum is not always 
very accurate especially for the high-spin part; for example, the calculated moment of inertia 
is almost constant or even decreases as a function of spin, while the observed inertia markedly 
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increases at high-spin states in most of the realistic cases. This is due to the fact that the 
rotational effect, e.g., the Coriolis anti-pairing effect, is not taken into account when only 
one intrinsic state is employed in the projection. However, if we intend to combine the 
cranking model and the angular momentum projection method, we encounter a puzzling 
situation; namely, we have many spectra depending on the cranking frequencies, each of 
which is obtained as a result of angular momentum projection from one mean-held state 
associated with each different frequency. It may be considered that the best way is to choose 
the lowest energy state for each spin value; i.e., the idea of the variation after projection 
regarding the cranking frequency as a variational parameter. It has been shown, however, 
that this does not necessarily work in practice [12|; if the rotational frequency is treated 
as a variational parameter at each spin value, it happens that the resultant spectrum does 
not compose a regular rotational pattern. In Ref. [12| it is recommended to calculate the 
spin I state by the projection from the cranked state with the frequency curot that gives 
(<hcr(<^rot)|-^i/|^cr(<^rot)) ~ / h. However, this procedure is not efficient because the angular 
momentum projection should be performed for each spin state from the intrinsic state with 
each different cranking frequency, and there is still no guarantee that a regular rotational 
sequence is obtained. 

Thus, how to compromise the cranking model and the angular momentum projection 
method is a non-trivial problem. In the present work, we attack the problem by a different 
approach: The resultant state for each spin value is calculated by superposing the multiple 
states obtained by projection from the cranked mean-held states with various rotational 
frequencies, i.e., in the sense of the conhguration mixing or the generator coordinate method 
(GCM). The original idea goes back to the work by Peierls-Thouless 13|. It will be shown 
that the resultant spectrum does not essentially depend on the chosen set of rotational 
frequencies as long as the number of them is sufficient; the necessary number is in fact 
rather small, something like four to hve. Thus the accurate rotational band can be efficiently 
calculated by the conhguration mixing, or the GCM, of a rather small number of angular 
momentum projected states; we call this procedure as “angular momentum projected multi- 
cranked conhguration mixing”. After a brief explanation of the formulation of the actual 
procedure in Sec. El we will show three examples of calculation in Sec. El Sec. |4] is devoted 
to summary and discussion. 


2 Multi-cranked configuration mixing 

The method of calculation we adopt is a version of the angular momentum projected 
conhguration mixing, or the projected GCM, where the cranking frequency Wrot is considered 
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as a generator coordinate. Thus the resultant state with spin I is obtained in the form 

= y^.9Kn,a ( 2 ) 

Kn 

where the operator is the angular momentum projector and the amplitude ^ is 

determined by the so-called Hill-Wheeler equation, 


'^Kn,K'n' 9K>n',a “ ^kn,K'n' 9k'n',ay ( 3 ) 

K'n' K'n' 

with dehnitions of the Hamiltonian and norm kernels, 


'^Kn,K'n' 

■^^Kn,K'n' 


($n 


H 

1 


KK' 




(4) 


see e.g. Ref. [2] for details. A set of the intrinsic states, 1$^), from which the projection 
is performed, is chosen to be composed of mean-held states, i.e., Hartree-Fock-Bogoliubov 
(HFB) states, which are obtained by the selfconsistent cranking model based on the Routhian 
H' in Eq. ([1]) with a suitably chosen set of rotational frequencies, ; n = 1, 2, • • • , Umax), 
i.e. 1$^) = |<hcr(i^roty • Then the state in Eq. ([2]) is nothing else but the discrete version of 


\'^M,a) — / hcUrot -PM_K-|4>cr(cUrot)), 


(5) 


K 


which was originally considered by Peierls-Thouless 13| . 

Instead of the amplitude gj^na hq. (Ej), the properly normalized amplitude j^. 


fk 


n,a 


9K'n',a-: 

K'n' 


( 6 ) 


should be considered in order to study the conhguration mixing. For example, the probability 
of the n-th HFB state in the eigenstate is given by 


PaW = (7) 

K 

or for Eq. ([5]), 

Pa(Wot) = l/i:,a(^rot)y (8) 

K 

In the present work we only consider the angular momentum projection for simplicity, 
but the number and parity projections should be incorporated if necessary. However, the 
conservation of neutron and/or proton number is important for the conhguration mixing 


4 



calculation, and we use the approximate way to take it into account [l^ when the system is 
in the superfluid phase; i.e., the Hamiltonian kernel in Eq. (I3]) is replaced by 


U 


Kn,K'n' 


H 


Kn,K'n' 


\y ANj^ 




(9) 


with 


A^L,A-n' = {^n\{N - = ($,|(Z - Zo)PkK'\^n'), (10) 

where the operator N (Z) is the neutron (proton) number operator and Nq (Zq) is the 
number to be flxed. The parameters Aiy and Ayr in Eq. ([9]) are taken to be the calculated 
neutron and proton chemical potentials of the first HFB state this choice is enough 

for the present purpose. 

In the following section we will show examples of calculation performed by the procedure 
of the angular momentum projected multi-cranked configuration mixing described above. We 
have recently developed an efficient method for the general quantum-number projection and 




GCM , and it has been applied to the study of the nuclear tetrahedral deformation [16| 


17l |. The method is fully utilized also in the present work. As for the effective interaction 


181 with the DIS parameter set 1^. The 


we employ the finite range Gogny interaction 
standard method of the harmonic-oscillator basis expansion is used for the HFB as well as 
the projection calculations. The calculational details related to the Gogny interaction are the 
same as those explained in Ref. 17|] except for one point: The Goulomb exchange contribution 
is treated in the Slater approximation although we are able to perform the exact calculation. 
This is because the selfconsistently calculated pairing correlation for protons quite often 
vanishes (or becomes very small) in the normal deformed nuclei in the rare earth region due 
to the fact that the Goulomb antipairing effect is too strong in the exact treatment 20|. This 


deficiency would be avoided if the variation after the number projection is performed 211. 


The proton pairing correlation is important to describe the collective rotation and we use 
the HFB calculation in the present work. Therefore the Slater approximation is the simplest 
way to avoid the deficiency. 


3 Examples of calculation 

In the present work we consider three examples of application; the ground state band of a 
typical rare-earth nucleus ^®^Er, the ground state band of a very neutron-rich nucleus ^®Mg, 
and the superdeformed band of a representative nucleus ^®^Dy. It has been found that the 
non-cranked HFB calculation gives axially symmetric deformation in these three examples. 
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We will discuss the A-pole deformation parameter of the calculated mean-held dehned by 




A 

4!r 


AR^ 

and the average pairing gap by j^| 

-J2^ab<b 


with R = 




' i=l 


A = 


a>b 


E 

a>b 


K 


ab 


with ^ah = E '^ab,cd b^cdi 

c>d 


( 11 ) 


( 12 ) 


where the quantity Kab is the abnormal density matrix (the pairing tensor) and is the 
matrix element of the pairing potential with the anti-symmetrized matrix element Vab,cd of 
the two-body interaction [^. 

As for the oscillator basis expansion the frequency huj = Al/A~^^^ MeV is employed 
and all the basis states with the oscillator quantum numbers {nx,ny,nz) satisfying Ux + 
Uy + Uz < are retained. We mainly use = 12; however, smaller values are also 

adopted for some calculations which require heavy numerical effort. For the calculation of the 
angular momentum projector in Eq. ([2]), the Gaussian quadrature over the three Euler angles 
(a,/ 5 , 7 ) is necessary: Even if the non-cranked HFB calculation gives the axially symmetric 
solution, the cranking term in Eq. ([T]) breaks the axial symmetry of the system. The number 
of integration mesh points for the angle a, which is the same for the angle 7 , is taken to 
be Na = Nj = 2JFinax + 2, and that for the angle P to be Np = 2Jniax + 2, where /max and 
//max are the maximum values of the I and K quantum numbers suitably chosen for each 


case. 


It has been known that the density-dependent term in the Skyrme interaction causes 


various numerical problems, see e.g. Refs. 23|, [2^ and references therein. The Gogny inter¬ 


action has the same density-dependent term and we encounter the situations that unphysical 
spectra are obtained due to the density-dependent term; for example, the spectra in the well- 
deformed odd mass nuclei sometimes do not look like rotational. The reason why we judge 
that it is due to the density-dependent term is that quite nice rotational spectra are recovered 
if we exclude it, although then the absolute energy and the moment of inertia are different 
from what they should be. It seems to us that the time-odd components of the wave func¬ 
tion induced by the cranking term tend to increase the rate of occurrence of this problem 
particularly for odd and odd-odd nuclei. However, we only consider even-even nuclei in the 
present work and the considered cranking frequency is not very high. Any problems related 
to the density-dependent term did not seem to occur in the following examples. 
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Another remark is related to the so-called norm cnt-off; the states whose norm eigenvalnes 
are smaller than a certain value are excluded when solving the Hill-Wheeler equation ([3]). 
The value for cut-off should be small enough not to miss important contributions, while it 
should not be too small in order to avoid the numerical difficulty in the GCM related to the 
vanishing norm states j^. We generally try to find a proper value as small as possible. It is, 
however, necessary to choose a suitable value for each calculation depending on the situation. 
This is because we use the same norm cut-off for all spin members under consideration in 
the calculation of a whole rotational band. It could happen that some small norm states 
included at some spins, which is controlled by the value of the norm cut-off, get excluded 
at different spins when solving Eq. (j3]). Such a change of included states at different spins 
may cause discontinuities of the rotational spectrum. Therefore we should adjust value of 
the norm cut-off so that any discontinuities would not happen; we are able to hnd such a 
proper value for the norm cut-off. 


3.1 Ground state band of 

The first example is the ground state rotational band of a nucleus ^®‘^Er, where the high- 
spin states up to = 22+ are measured 2^. We use the oscillator basis with = 12. 

The selfconsistently calculated values of the mean-field parameters of the ground state in this 
nucleus are (32 = 0.311 for the quadrupole deformation, which roughly corresponds to the 
experimentally deduced value 2^, and A = 0.874, 0.906 MeV for the neutron and proton 
average pairing gaps, respectively, which are about 10%-13% smaller than the even-odd mass 
differences. In this case the system is in the superfluid phase, and we apply the quasiparticle 
basis cut-off at the level of 10“®, i.e., only the canonical basis states whose occupation 
probabilities are larger than 10“® are considered, see Ref. 1^ for details. As for the projection 


calculation the maximum values of the I and K quantum numbers are taken to be /max = 22 
and //max = 16 and therefore the numbers of integration mesh points are Na = = 34 and 

Np = 46 for the Euler angles. 

First we show the calculated excitation spectrum obtained by the angular momentum 
projection from the non-cranked HFB state (i.e. curot = 0) in comparison with the experimen¬ 
tal data in Fig. [T] The energy gain of the ground state by the angular momentum projection 
is 2.97 MeV. The value of norm cut-off can be taken to be 10“^^ in this calculation. As it is 
clearly seen in the figure, the projected rotational energies from the non-cranked state are 
systematically larger than the experimental ones. We also include the result of projection 
from the cranked HFB state with a small cranking frequency hw^ot = 0.01 MeV in Fig. [T] The 
curvature of the spectrum as a function of spin decreases considerably as a result of including 
the time-odd components in the wave function by the cranking model, and the agreement 


7 









Fig. 1 Excitation energy spectra of the ground state rotational band in ^®'^Er calculated 
by the projection method from the non-cranked HFB state (the label “curot = 0.0”) and from 
the cranked HFB state with hcurot = 0.01 MeV (“curot = 0.01”). Experimentally measured 
spectrum (“Exp.”) is also included. The O"*" energy is taken as the energy origin in each 
spectrum. The oscillator basis with = 12 is used. 


with the experimental data is much better. Namely the calculated moment of inertia is 
enlarged by the small cranking term nearly by a factor of two. This result clearly shows 
the importance of the cranking procedure for the angular momentum projection calculation, 


which was emphasized in Refs. 1^, ll6| . However, the rotational energies are overestimated 


at the higher spin states, / > 16, although they are slightly underestimated at the lower spin 
states, I < 12: The experimental moment of inertia increases at higher-spin states, while the 
calculated inertias are almost constant as long as the projection is performed from only one 
intrinsic HFB state, see also Figs. Eland[3 Thus it is necessary to improve the calculation 
by the multi-cranked configuration mixing. 


We show in Fig. |2]the results of the projected multi-cranked conhguration mixing and of 
the simple projection from only one cranked HFB state with hcurot = 0.01 MeV in compari¬ 
son with the experimental data. The reference rotational energy /(/ + l)h^/(2j7’) with JT" = 
40 [h^/MeV] is subtracted to display the detailed behaviors more clearly. A set of five almost 
equidistant cranking frequencies, n = 1 ; 5) = (0.01, 0.05, 0.10, 0.15, 0.20) MeV, is 

adopted. The agreement with the experimental data is much better; the difference of exci¬ 
tation energy is less than 130 keV in the whole spin range, 0 < / < 22, if the configuration 
mixing is performed. In order to display the effect of the configuration mixing, the five spec¬ 
tra obtained by the projection from one intrinsic HFB state with the hve different cranking 
frequencies are shown in Fig. [3] in addition to the result of the projected configuration mix¬ 
ing. The energy origin is chosen to be the O'*" energy of the final configuration mixing and 







Fig. 2 Excitation energy spectra subtracting the reference rotational energy, /(/ + 
l)/80 MeV, for ^®^Er. The result of the simple projection from one cranked HFB state 
with hjjJrot = 0.01 MeV as well as that of the projected multi-cranked conhguration mixing 
(the label “Mixed”) explained in the text are included in addition to the experimental data. 
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Fig. 3 Energy spectra for ^®'^Er obtained by the simple projection from one intrinsic HFB 
state with Eve values of the cranking frequencies, hoj^ot = 0.01, 0.05, 0.10, 0.15, 0.20 MeV, 
compared with the result of the projected configuration mixing employing those hve HFB 
states. The energy origin is taken as the energy of the 1 = 0 state of the conhguration mixing 
calculation, and the same reference rotational energy is subtracted as in Fig. [2l 


the reference rotational energy is subtracted as in the same way as in Fig. O The resultant 
rotational spectra are similar as long as one cranked HFB state is used, while the absolute 
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O"*" energies are smaller when projected from the HFB state with higher cranking frequencies. 
The energy gain of the conhguration mixing from the simple projection with the non-cranked 
HFB state is about 0.83 MeV and the total energy gain by the projected conhguration mix¬ 
ing from the HFB energy is 3.81 MeV. It is quite interesting to mention that the excited 
K'^ = 0"*" rotational band of the Hill-Wheeler equation as a result of the conhguration mix¬ 
ing has about 6 MeV higher energy than the ground state. This is consistent to the fact 
that there is only one rotational band associated with the ground state, and the hve rota¬ 
tional sequences obtained by the simple projection from each conhguration are not totally 
independent; in fact the overlaps between these hve HFB states are rather large. 



Fig. 4 Probability distribution (Eq. ([H])) over the hve HFB conhgurations with diherent 
cranking frequencies for the spin I member of the ground state rotational band of ^®'^Er 
obtained by the projected conhguration mixing. 


The probabilities dehned in Eq. ([8]) of the hve conhgurations associated with hoj^ot = 
0.01, 0.05, 0.10, 0.15, 0.20 MeV are depicted in Fig. 01 The distributions shown in Fig. 0]are 
rather diherent from what are expected, i.e., the probabilities have peaks at the cranking 
frequencies corresponding to (4’cr(i^rot)|>^i/|*^’cr(Wot)) ~ / h. In fact the distributions spread 
over hve conhgurations with the middle one (hcurot = 0.10 MeV) having always the lowest 
probability. In fact the distributions for all the 1 = 0,4, 8, 12 members are rather similar and 
the main diherences are the probabilities of the hrst (hajrot = 0.01 MeV) and the last (hwrot = 
0.20 MeV) conhgurations. It should be pointed out that the rotational sequence projected 
from the last conhguration is lowest in energy, see Fig. [3l and therefore its probabilities are 
relatively large for all the / = 0, 4, 8, 12 members, which may be specihc in this example. 

Figure [5] depicts the angular momenta as functions of the cranking frequency for various 
calculations in comparison with the experimental data. For the cranked HFB calculation the 
angular momentum is calculated by /(wrot) ^ = (‘^’cr(Wot)|'^y|‘^’cr(Wot)) — as a function 
of the rotational frequency curot, where the subtraction of is the semiclassical correction. 
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Fig. 5 Angular momentum versus rotational frequency for ^®'^Er obtained by various 
calculations; the cranked HFB (the label “CHFB”), the simple projection from the non- 
cranked HFB state (“curot = 0.0”), and the projected conhguration mixing (“Mixed”) in 
comparison with the experimental data (“Exp.”). 


see e.g. Ref. Q. For the results of projection calculation either from the non-cranked HFB 
state or with the configuration mixing, the rotational frequency is calculated by the finite 
difference, hjJxot{I) = {,E{I + 1) — E{I — l))/2, where E{I) is the energy at spin I. One 
can see that the slopes at lower spins are quite different for various calculations and again 
the agreement of the configuration mixing calculation with the experimental data is best, 
while the slope of the cranked HFB is larger and that of the simple projection from the 
non-cranked HFB state (hcurot = 0.0) is considerably smaller than the experimental data. 
It should be noticed that the abrupt upbend of the angular momentum is observed in the 
cranked HFB calculation at hcurot ~ 0.25 MeV. This is caused by the angular-momentum- 
alignment of the lowest two quasi-neutrons along the rotational axis, which leads to the 
crossing between the ground state band (g-band) and the Stockholm band (s-band). This 
band-crossing has been known for many years as the o^in of the backbending phenomenon 
hrst observed in Ref. 27|, see e.g. discussions in Refs. [71, l 8 |. The effect of this band-crossing 
is reflected as an irregularity at / = 15 in the experimental data in Fig. O while there is no 
such effect for the projected calculations because the crossing with the s-band configuration 
is not considered. 


It should be mentioned that the selfconsistent cranked HFB calculation shows the 
characteristic changes of the mean-field parameters as functions of the rotational fre¬ 
quency. For example, the quadrupole deformation parameter ^2 slightly increases until about 
hcjrot ~ 0.2 MeV and starts to decreases, although the total amount of change is less than 10% 
in the depicted range of frequency in Fig. [5l The triaxiality parameter 7 gradually increases 
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from zero at zero frequency up to a small value about 7 +7° at hw^ot ~ 0.24 MeV. The 

average pairing gaps for both neutrons and protons gradually decrease and that of neutrons 
suddenly drops after the g-s crossing because of the blocking effect. These behaviors of the 
mean-held parameters have been well-known for well deformed nuclei in the rare-earth region 
by similar selfconsistent cranking calculation but with simple schematic interactions, see e.g. 
Ref. for this nucleus. In fact the results of the present cranked HFB calculation are very 


similar to those of Ref. 


(note that the cranking axis is the |/-axis in the present work 


while it is the x-axis in Ref. |28|]) except for the fact that the neutron pairing gap vanishes 
after the crossing frequency fvjJrot ~ 0.24 MeV in the present calculation. This is the main 
reason why we do not investigate the s-band in the present work: The calculated moment of 
inertia for the s-band is overestimated due to the vanishing neutron pairing correlation. We 
have used the maximum cranking frequency, /kUrot = 0.20 MeV, for the conhguration mixing 
calculation of the g-band because of a similar reason; the effects of the s-band conhguration 
should be excluded for describing the g-band. 

The present cranked HFB calculation in Fig. |5]does not show the backbending behavior 


like in Ref. 


, because the minimum-search with the proper treatment of the angular 


momentum constraint [2j is not performed in the backbending region. The proper treatment 
of the band-crossing phenomenon within the projected conhguration mixing approach is now 
under investigation. 



Fig. 6 Moments of inertia versus spin value for ^®'^Er obtained by various calculations in 
comparison with the experimental data like in Fig. |5l see the text for the precise dehnition 
of the hrst moment of inertia 
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Fig. 7 Moments of inertia versus spin value for ^®'^Er obtained by the simple pro¬ 
jection from one intrinsic HFB state with hve values of the cranking frequency, huJxot = 
0.01,0.05,0.10,0.15,0.20 MeV, compared with the result of the projected conhguration 
mixing employing those hve HFB states. The ordinate scale is enlarged from that of Fig. [Ul 


In order to study the property of the rotational band more closely, we depict the 
result of the so-called hrst (or kinematic) moment of inertia for various calcula¬ 

tions and the experimental data in Figs. |6] and [71 It is dehned for the cranked HFB 
calculation by = (‘hcr(i^rot)|Ty|<hcr(i^rot))/i^rot, which is plotted as a function of 

I fi = (*hcr(i^rot)|'7y|*hcr(ci;rot)) “ and for the projection calculations and the experimen¬ 
tal data by = {21 + l)h^/(F(/ -1-1) — E{I — 1)). As it is clearly seen in Fig. El the 

calculated moment of inertia obtained by the simple projection from the non-cranked HFB 
state is considerably underestimated: Its value is only 66% of the experimental value at the 
lowest spin, while the results of the cranked HFB and of the simple projection from one 
cranked HFB state with the hnite frequencies hwxot < 0.20 MeV are about 15% larger than 
the experimental value at the lowest spin, which may be due to the fact that the calcu¬ 
lated pairing gaps are about 10% smaller than the even-odd mass differences. Compared to 
this, the result of the projected conhguration mixing gives slightly better agreement with the 
experimental data. It is interesting to notice that all the results of the simple projection from 
one cranked HFB state associated with different frequencies are rather similar at the lowest 
spin, see Fig. [71 As for the spin-dependence of the moment of inertia, only the results of the 
cranked HFB and of the projected conhguration mixing increase as functions of spin simi¬ 
larly to the experimental moment of inertia. The increase of calculated inertia by the cranked 
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HFB method is mainly caused by the gradual reduction of pairing correlations, i.e., the so- 
called Coriolis anti-pairing effect. Those of the simple projection from one HFB state with 
various cranking frequencies (including zero-frequency) are rather constant, or even slightly 
decreasing for the cases of the hnite frequencies, huj^ot = 0.01,0.15,0.10,0.15 MeV, and are 
completely different from the experimental spin-dependence. Thus the configuration mixing 
is important in order for the projected result of moment of inertia to increase, although the 
amount of increase is not enough compared with the behavior of the experimental data. 

We have used a set of almost equidistant cranking frequencies to generate the cranked 
HFB states for the configuration mixing. There are, however, inhnitely many possible choices. 
It is interesting to see how the result of conhguration mixing depends on the chosen set of 
cranking frequencies, which is analyzed in Fig. [H] by making use of random numbers. In 
this calculation the smaller basis size with = 10 is used to reduce the numerical task. 

We have randomly generated various sets of cranking frequencies with the conditions that 
the frequency satisfies, 0 < hwrot < 0.20 MeV, to avoid the effect of the g-s band crossing, 
and the difference of two nearest frequencies satisfies, hAcUrot > 0.02 MeV, to avoid too 
large overlap between the two associated HFB states, which may cause the vanishing norm 
problem of the GCM procedure j^. It happens that the norm cut-off problem occurs for 
these randomly generated frequencies. Therefore, we had to adjust values of norm cut-off in 
the range 10“® — 10“^®. The results of analysis employing three, four, and hve conhgurations 
associated with randomly chosen frequencies are shown in the first, second, and third rows, 
respectively, in Fig. [HI We have performed eight trial calculations in each case, and the average 
and standard deviation of the calculated moments of inertia are plotted in the left panels, 
where the standard deviation is shown as error bars. The middle panels show all moments 
of inertia calculated in eight trials, and the randomly generated sets of cranking frequencies 
applied in these trials are displayed in the right panels. As is seen in Fig. [HI when the number 
of conhgurations is three, not only the standard deviation is large but also the average value 
itself deviates considerably from the one calculated with larger numbers of conhgurations. 
The diherence of calculated moments of inertia among eight trials is getting smaller as the 
number of mixed conhgurations is increased, even though the generated sets of frequencies 
are markedly diherent. In fact, if the hve conhgurations are employed the standard deviations 
are within the size of symbols, see the left-bottom panel in Fig. [HI Therefore the result of 
multi-cranked conhguration mixing with the equidistant set of hve cranking frequencies gives 
the almost unique result. The small diherence between the converged result of mixing the 
hve conhgurations in Fig. [H]and that in Fig. [6] or Fig. [7|is due to the diherent model space 
with = 10 and 12 being used. 
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Fig. 8 Analysis of moments of inertia for ^®'^Er calculated by the projected configura¬ 
tion mixing employing cranked HFB states with randomly chosen cranking frequencies. The 
smaller model space = 10 is used in this calculation. Three sets of three panels in 

the hrst, second, and third raw show the results of calculation with mixing three, four, and 
hve conhgurations, respectively, in comparison with the experimental data. In each set the 
left panel depicts the average moment of inertia with the error bars indicating the standard 
deviation for eight trials of randomly chosen frequencies, which are shown in the right panel 
as crosses, and the middle panel displays the moments of inertia calculated in those eight 
trials. The results of eight trials are distinguished by different symbols. 


At the end of this subsection we briefly discuss the transition probabilities. We compare 
the calculated in-band B{E2) values with experimental data [2^ in Fig. [HI a good agreement 
is achieved. In order to show the effect of conhguration mixing for the B{E2) values, we 
show those calculated by the simple projection from one HFB state with different cranking 
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Fig. 9 Calculated in-band B{E2; 1^1 — 2) values by the projected configuration mixing 
in for the ground state band of ^®‘^Er in comparison with the experimental data. 



Fig. 10 Calculated in-band B{E2; 1^1 — 2) values by the simple projection from one 
HFB state with different cranking frequencies in comparison with the result of configuration 
mixing. The ordinate scale is enlarged compared with Fig. [9l 


frequencies in Fig. [101 The calculated B{E2) value is larger when projected from the HFB 
state with higher frequency; this can be well understood because the deformation parameter 
P 2 slightly increases and the triaxiality parameter 7 is positive and also slightly increases in 
this frequency range (0 < < 0.2 MeV), both of which lead to the increase of quadrupole 

moment around the rotation axis (note that it is the y-axis in the present work). As it is 
shown in Fig. [101 however, the configuration mixing makes the B{E2) values smaller than any 
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results of the simple projection from one HFB state at all spin values considered, although 
the amount of reduction is rather small, less than 5%. We think that this is due to the fact 
that the mixing probabilities rather spread over hve conhgurations as was shown in Fig. lU 
which reduces the expectation value of the E2 operator by decoherence. As it is well known 
that the B[E2) in the rotational band can be nicely htted by the rotor model [l|: In fact the 
quadrupole moment is usually extracted from the experimental data assuming this model. 
The angular momentum projection calculation of the B{E2) for well-deformed nuclei justihes 
the rotor model, see e.g. Ref. j^. Since the deformation parameter of the selfconsistent HFB 
calculation almost reproduces the correct quadrupole deformation, the agreement in Fig. |9] 
may be natural. 


3.2 Ground state band of 

We take a nucleus as an example of the deformed unstable nuclei; It has been 

predicted, see e.g. Ref. j^, that the N = 28 Mg isotope is deformed in spite of the fact 
that = 28 is one of the magic numbers in stable nuclei. This nucleus is also predicted to 
be the drip-line nucleus in Ref. 30|. We use the oscillator basis with = 12 for this 

example. The selfconsistently calculated value of the quadrupole deformation parameter is 
(32 = 0.334, and the pairing correlations for both neutrons and protons vanish. The neutron 
skin develops considerably in this nucleus and the skin thickness is near 1 fm along the 
short-axis of deformation. The calculated root-mean-square radius is 3.63 fm, which is about 
15% larger than the value given by the empirical radius 1.2A^/^ fm. It should be noted 
that the use of the calculated radius R in the dehnition fITTll is important to estimate the 
deformation; the deformation parameter is largely overestimated if the value R = 1.2A^/‘^ fm 
is employed instead. These basic features of the present result of the Gogny HFB calculation 
well correspond to those of Ref. j^|, although the calculated deformation is slightly larger 
in our calculation. As for the projection calculation we take /max = 20 and //max = 12 and 
therefore the numbers of integration mesh points are = Aly = 26 and Np = 42 for the 
Euler angles. The preliminary report of our calculation for this nucleus was published in 
Ref. [31]. 


The excitation spectra calculated by the angular momentum projection method are 
shown in Fig. [HI where are shown three spectra obtained by the simple projection from 
the non-cranked HFB state and from the cranked HFB state with a small frequency 
hcurot = 0.01 MeV, and by the projected conhguration mixing. A set of four equidistant 
cranking frequencies, n = 1 : 4) =(0.01, 0.34, 0.67, 1.00) MeV, is adopted for calcu¬ 

lating the cranked HFB states employed in the conhguration mixing calculation. There is no 
experimental data available yet for this nucleus. The value of norm cut-off is taken to be 10“® 
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Fig. 11 Excitation energy spectrum of the ground state rotational band in ^*^Mg calculated 
by the projected configuration mixing (the label “Mixed”) in comparison with those by the 
simple projection from the non-cranked HFB state (“curot = 0.0”) and from the cranked HFB 
state with hu^ot = 0.01 MeV (“curot = 0.01”). The oscillator basis with = 12 is used. 
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in this calculation. The energy gain of the ground state by the simple angular momentum pro¬ 
jection is 3.07 MeV and that by the configuration mixing is 3.47 MeV. The projected GCM 
calculation for this nucleus with the same Gogny DIS interaction was reported in Ref. j^. 
Although the detailed treatment is different from the present calculation, the results are 
quite consistent with each other. 

The effect of the conhguration mixing is shown in Fig. [121 where the four spectra 
obtained by the projection from one intrinsic HFB state with four cranking frequencies, 
l^rot = 0.01, 0.34, 0.67, 1.00 MeV, are depicted in addition to the result of the projected con¬ 
figuration mixing employing those cranked HFB states. The probabilities defined in Eq. ([H]) 
of these four configurations are also shown in Fig. [T31 The energy origin is chosen to be the O'*" 
energy of the final configuration mixing and the reference rotational energy, I{I + 1)/14 MeV, 
is subtracted. The resultant rotational spectra of the projection from one cranked HFB state 
are similar at lower spins but considerably different at higher spins. The absolute O”*" energies 
are larger when is used the HFB state with higher cranking frequency, which is opposite to 
the case of the previous example ^®‘^Er. The excited K'^ = 0"*" rotational band in the conhgu¬ 
ration mixing has about 2.7 MeV higher energy than the ground state, and again it indicates 
that there is only one rotational band associated with the ground state. The probability dis¬ 
tributions in Fig.Hniseem to be consistent with what are expected, i.e., the probabilities have 
peaks roughly at the cranking frequencies corresponding to ($cr(<^rot)|<7y|4>cr(i^rot)) ~ 7 h., see 
also the following Fig. [141 
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Fig. 12 Energy spectra for obtained by the simple projection from one intrinsic 

HFB state with fonr valnes of the cranking freqnency, hwrot = 0.01, 0.34, 0.67, 1.00 MeV, 
compared with the resnlt of the projected conhgnration mixing employing those fonr HFB 
states. The energy origin is taken as the energy of the 1 = 0 state of the conhgnration mixing 
calculation, and the reference rotational energy, /(/ + 1)/14 MeV, is subtracted. 



TlWrot [MeV] 


Fig. 13 Probability distribution (Eq. (0)) over the four HFB conhgurations with different 
cranking frequencies for the spin I member of the ground state rotational band of ^'^Mg 
obtained by the projected conhgnration mixing. 


Figured!] depicts the angular momenta as functions of the cranking frequency for various 
calculations. One can see that the slope of the result of the projection from the non-cranked 
HFB state is smaller than the others, while the three results of the cranked HFB, of the 
simple projection from the cranked HFB state with hw^ot = 0.01 MeV, and of the projected 
conhgnration mixing are rather similar at lower frequencies. However, at higher frequencies, 
hcjrot ~ 1.3 MeV, the upbending behavior can be seen for the projected conhgnration mixing 
calculation. We do not understand the precise reason but it seems that the rotational band at 
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Fig. 14 Angular momentum versus rotational frequency for obtained by various 

calculations; the cranked HFB (the label “CHFB”), the projection from the non-cranked 
HFB state (“curot = 0.0”), the projection from the cranked HFB state with = 0.01 MeV 
(“curot = 0.01”), and the projected conhguration mixing (“Mixed”). 


higher spins, / > 12, changes its nature. This behavior does not disappear even if the cranked 
HFB states with higher frequencies, hcjrot >1.0 MeV, are included in the configuration 
mixing calculation. Taking a closer look at the result of the cranked HFB in Fig. [HI the 
slope gradually increases at the higher frequency, hw^ox ^1-2 MeV, which is caused by the 
gradual rotational alignment. The last occupied neutron orbit [321 is bound but its binding 
energy, 1 MeV at Wrot = 0, decreases to about 150 keV at higher frequency hcurot ~ 1.4 MeV, 
where this orbit strongly interacts with a discretized continuum state; note that the pairing 
correlations vanish in this nucleus. Therefore we only consider the rotational spectrum in 
the spin range, I < 12. 



Fig. 15 Moments of inertia versus spin value for ^*^Mg obtained by various calculations 
like in Fig. dH 
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Fig. 16 Moments of inertia versus spin value for obtained by the simple pro¬ 

jection from one intrinsic HFB state with four values of the cranking frequency, hcurot = 
0.01,0.34,0.67,1.00 MeV, compared with the result of the projected conhguration mixing 
employing those four HFB states. 


In Figs. [15] and dB] we compare the results of the moment of inertia obtained by 
various calculations as functions of spin. The inertia by the simple projection from the non- 
cranked HFB state is smaller than the others and clearly shows the importance of the effect 
of cranking to increase the moment of inertia. However, the effect is not so conspicuous as 
in the case of ^®'^Er; the calculated moment of inertia by the simple projection from the 
non-cranked HFB state is about 80% of that by the projected conhguration mixing. It is 
interesting to mention that all inertias calculated by the simple projection from one HFB 
state decrease as functions of spin in this nucleus; see Fig. dBl Especially the reduction of 
the result with the smallest cranking frequency, fvjJrot = 0.01 MeV, is largest. This kind of 
behavior is quite unusual and it was speculated in Ref. that this reduction of inertia 
may be characteristic for the unstable skin nucleus. As it is shown in Figs. [I5] and fTOl 
however, the conhguration mixing recovers the usual behavior of the moment of inertia: It 
is almost constant or gradually increases at higher spins. We believe that the result of the 
conhguration mixing is more reliable and the decrease of the moment of inertia as a function 
of spin 3l| would not be realized. Although the reduction of inertia seems to be unrealistic, 
it is still interesting to investigate the rotational motion in the deformed unstable nuclei up 
to the neutron drip-line. A systematic Gogny HFB and the projection calculations for the 
Mg isotopes including '^'^Mg are now in progress [^ . 
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Fig. 17 Analysis of moments of inertia for calculated by the projected configura¬ 

tion mixing employing cranked HFB states with randomly chosen cranking frequencies. The 
smaller model space = 8 is used in this calculation. Two sets of three panels in the 

hrst and second raw show the results of calculation with mixing two and four conhgurations, 
respectively. In each set the left panel depicts the average moment of inertia with the error 
bars indicating the standard deviation for eight trials of randomly chosen frequencies, which 
are shown in the right panel as crosses, and the middle panel displays the moments of iner¬ 
tia calculated in those eight trials. The results of eight trials are distinguished by different 
symbols. 


In order to see how the result of conhguration mixing depends on the choice of cranking 
frequencies, similar analysis to the case of ^®'^Er has been performed for ^*^Mg by using 
randomly chosen frequencies; the result is shown in Fig. [T71 In this calculation the smaller 
basis size with = 8 is used to reduce the numerical task because the conhguration 

mixing calculations should be repeated many times. The random numbers in this case are 
generated in the range, 0 < hcurot <1-0 MeV, and with the condition that the difference of 
two nearest frequencies satishes, hAcUrot >0.1 MeV, to avoid the large overlap between the 
two associated HFB states. The results with the two and four conhgurations are shown in 
the hrst and second raw, respectively, in Fig. [T71 The eight sets of cranking frequencies are 
generated for trials, which are shown in the right panels of the hgure, and the eight moments 
of inertia calculated by the projected conhguration mixing in each trial are displayed in the 
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middle panels. The average of these eight trials is shown in the left panels with the error 
bars being the standard deviation; in fact the standard deviations are rather small except 
at the highest spin. As it is clear, the mixing with four configurations gives rather similar 
results for all the eight trials, which is in contrast to the case of ^®'^Er, where the standard 
deviations are small only in the result with mixing five configurations. The difference between 
the converged result of mixing the four conhgurations in Fig. [T3 and that in Fig. [T5]is due to 
the different model space with = 8 and 12 being used. The large standard deviations 

at the highest spin in this example may be related to the irregularity corresponding to the 
gradual alignments at higher rotational frequency in Fig. [141 


3.3 Superdeformed band of 

The last example is the superdeformed band of a nucleus ^®^Dy. It is the first nucleus, 
in which the high-spin superdeformed rotational band was discovered by the discrete-line 
gamma-ray spectroscopy 3^. Moreover, the linking transitions between the superdeformed 


and normal deformed states were measured in this nucleus 


35l | , and the excitation energy and 


the precise spin-assignment of the superdeformed band are known. We have done the HFB 
and the projection calculations using the oscillator basis with = 12 for this example. It 

has been found that the superdeformed minimum exists in the non-cranked HFB calculation, 
where the neutron pairing vanishes but the very weak proton pairing remains with the 
average pairing gap ~ 0.4 MeV. The projection calculation with the pairing correlation 
is much more time consuming and we need to calculate up to very high-spin states like 
I ^ 60. Therefore we switched off the proton pairing correlation for the calculation of this 
example. The difference between the minimum of the binding energies with and without 
the proton pairing correlation is less than 20 keV within the mean-field approximation. The 
selfconsistent deformation is axially-symmetric and the quadrupole deformation parameter is 
132 = 0.715, which decreases gradually up to 0.699 when the system is cranked up to = 
0.7 MeV keeping very well the axial-symmetry. This result of the mean-field calculation 
well corresponds to those of the Woods-Saxon- or Nilsson-Strutinsky calculations of, e.g.. 
Refs. 


36 


37l |. although the magnitude of deformation is slightly larger than the results of 


these references. This difference is due to the different definition of the deformation; in the 
Strutinsky calculation the deformation is defined for the avera ge p otential, while in the 
HFB calculation it is defined by the density distribution, see Ref. (36| . In fact the calculated 
deformation well reproduces the measured quadrupole moment, see the end of this subsection 
related to Fig. [211 As for the projection calculation we take /max = 62 and //max = 22 and 
therefore the numbers of integration mesh points are = N.y = 46 and Np = 126 for the 
Euler angles. 
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In this example we do not show either the rotational spectrum nor the angular momentum 
versus rotational frequency relation, but rather we concentrate on the moment of inertia. The 
projected configuration mixing calculation has been performed with four cranked HFB states 
with the equidistant cranking frequencies, n = 1 : 4) =(0.01, 0.24, 0.47, 0.70) MeV. 

with the value of norm cut-off taken to be 10“®. The energy gain of the 0"*" superdeformed 
state by the simple projection from the non-cranked HFB state is 4.11 MeV and that by the 
projected conhguration mixing is 4.14 MeV; the difference is rather small in this case. The 
calculated excitation energy of the superdeformed 24”'' state by the configuration mixing is 
9.97 MeV, which is compared to the experimentally measured value 10.644 MeV SSj; the 
agreement is rather nice. 



Fig. 18 Moments of inertia versus spin value for ^^^Dy obtained by the two calculations, 
the projection from the non-cranked HFB state (“wrot = 0.0”) and and the projected config¬ 
uration mixing (“Mixed”) in comparison with the experimental data (“Exp.”). The oscillator 
basis with N^ax = 12 is used. 


The moments of inertia for various calculations and experimental data are displayed in 
Figs. [IS and UHl separately. As is clearly seen, all the calculated inertias are rather constant 
as functions of spin in accordance with the experimental data, although the absolute value 
of the result of the projected configuration mixing is slightly larger than the experimental 
one. On the other hand the result of the simple projection from the non-cranked HFB state 
is considerably smaller than the experimental one, and again indicates the importance of the 
time-odd components in the wave function induced by the cranking term. The semiclassical 
moment of inertia of the cranked HFB calculated in the same way as in the previous examples 
coincides very well with that of the projected configuration mixing in this example; the 
ordinate scale is extremely enlarged to see the difference in Fig. [191 
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Fig. 19 Moments of inertia versus spin value for obtained by the cranked HFB (the 

label “CHFB”) is compared with the result of the projected configuration mixing (“Mixed”). 
The ordinate scale is enlarged compared with Fig. [THlto see the difference. 



Fig. 20 Moments of inertia versus spin value for ^®^Dy obtained by the simple 
projection from one intrinsic HFB state with four values of the cranking frequencies, 
hcjrot = 0.01, 0.24, 0.47, 0.70 MeV, respectively, compared with the result of the projected 
conhguration mixing employing those four HFB states. 


In order to show the effect of configuration mixing, we depict in Fig. [20] four moments 
of inertia obtained by the simple projection from one cranked HFB state with different 
frequencies, hujxot = 0.01, 0.24, 0.47, 0.70 MeV, in addition to the result of the projected 
conhguration mixing. Interestingly enough the calculated inertias by the simple projection 
are considerably different, both their absolute values and spin-dependences, as is clearly seen 
from the hgure; the one associated with the higher cranking frequency has larger absolute 
value and more strongly decreases as a function of spin. This is in contrast to the previous 
two examples of the ground state bands in ^®'^Er and ^*^Mg. Although the strong reduction of 
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Fig. 21 Probability distribution (Eq. ([H])) over the four HFB configurations with different 
cranking frequencies for the spin I member of the superdeformed rotational band of 
obtained by the projected conhguration mixing. 


these inertias completely contradicts with the experimental data, the projected conhguration 
mixing makes the resultant inertia almost constant in accordance with the tendency of 
the experimentally measured inertia. This fact clearly shows again the importance of the 
conhguration mixing to obtain the correct rotational behavior of the superdeformed states 
by the angular momentum projection method. In order to study how the four conhgurations 
mix, the probability distributions are depicted in Fig. [2T] In this case the distributions are 
just what are expected, i.e., the probabilities have nice peaks at the cranking frequencies 
corresponding to (‘hcr(i^rot)|</i/|d>cr(i^rot)) ~ / h, in contrast to the case of ^®^Er. This clearly 
shows that the superdeformed rotational band can be interpreted as a very good semiclassical 
rotor. 
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Fig. 22 The second moment of inertia versus spin value for ^®^Dy obtained by the pro¬ 
jected conhguration mixing (“Mixed”) in comparison with the experimental data (“Exp.”); 
see the text for the precise dehnition of 
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Figure |22] depicts t 
moment of inertia 


le result of configuration mixing for the so-called second (or dynamic) 
defined by = Ah^/{E{I + 2) + E{I - 2) - 2E(/)), in com¬ 

parison with the experimental data. The spin-assignment is not necessary to estimate this 
inertia and therefore it has been used frequently for the study of high-spin superdeformed 
bands. As it is seen in the hgure, a better agreement with the experimental data is obtained 
compared with the hrst moment of inertia in Fig. [181 although the calculated inertia is still 
slightly overestimated. 
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Fig. 23 Analysis of moments of inertia for ^®^Dy calculated by the projected conhgura- 
tion mixing employing cranked HFB states with randomly chosen cranking frequencies. The 
smaller model space = 10 is used in this calculation. Two sets of three panels in the 

hrst and second raw show the results of calculation with mixing two and four conhgurations, 
respectively, in comparison with the experimental data. In each set the left panel depicts the 
average moment of inertia with the error bars indicating the standard deviation for eight 
trials of randomly chosen frequencies, which are shown in the right panel as crosses, and the 
middle panel displays the moments of inertia calculated in those eight trials. The results of 
eight trials are distinguished by different symbols. 


It is interesting to see how the result of conhguration mixing depends on the choice 
of cranking frequencies also in this case. Similar analysis to the cases of ^®'^Er and ^*^Mg 
by employing the randomly chosen sets of cranking frequencies has been performed and 
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the result is shown in Fig. [231 In this calculation the smaller basis size with = 10 

is used to reduce the numerical task. The random numbers in this case are generated in 
the range, 0 < hwrot <0.7 MeV, and with the condition that the difference of two nearest 
frequencies satisfies, hAwrot > 0.07 MeV. The eight sets of cranking frequencies are employed, 
and the results with the two and four configurations are shown in the first and second raw, 
respectively, in Fig. |23l The left panels depict the average of eight trials with the error bars 
being the standard deviation, the middle panels the eight moments of inertia calculated in 
each trial, and the right panels the generated sets of frequencies. Surprisingly the mixing 
with only two conhgurations already gives a rather converging result, and the use of four 
configurations is enough to attain the almost completely unique result: Note that the ordinate 
scale is very enlarged in Fig. (231 compared with Fig. [THl The error bars are invisible in the left- 
bottom panel and the results of eight trials cannot be distinguished even though the eight sets 
of randomly chosen frequencies are considerably different. This result clearly conhrms again 
that a rather small number of configurations, four in this case, is sufficient to obtain reliable 
results. The difference of about 2% between the result of mixing the four configurations in 
Fig. 123] and that in Fig. [I8]is due to the different model space with = 10 and 12 being 
used. 
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Fig. 24 Calculated in-band B{E2\ 1^1 — 2) values by the projected conhguration mix¬ 
ing in for the superdeformed state band of ^®^Dy in comparison with the experimental 
data. 
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Finally we compared the calculated B{E2) values with the experimental data [39| in 
Fig. 1231 Although the increasing trend of the experimental data at higher spins, I > 50, is not 
reproduced, nice agreement is obtained. This indicates that the selfconsistently calculated 
quadrupole deformation agrees very well with the experimental data at least at the spin 
region, I < 50. In order to show the effect of conhguration mixing for the B{E2) values, those 




Fig. 25 Calculated in-band B{E2] / — 2) values by the simple projection from one 

HFB state with different cranking frequencies in comparison with the result of conhguration 
mixing. The ordinate scale is enlarged compared with Fig. [241 


calculated by the simple projection from one HFB state with different cranking frequencies 
are depicted in Fig. [25l As it clear from Fig. [25l the slight decrease of the calculated B{E2) 
in the spin range, I > 30, is a result of the conhguration mixing. This can be understood as 
follows: The selfconsistently calculated deformation parameter 132 decreases as a function of 
the cranking frequency with keeping very well the axial-symmetry as was already mentioned, 
which leads to the result that the B{E2) value calculated by the simple projection is smaller 
when projected from the HFB state with higher frequency as is also shown in Fig. [23 The 
mixing probabilities in this case have sharp peaks at the frequencies corresponding to the 
spin values under consideration, see Fig. [211 Therefore the quadrupole moment calculated 
by the conhguration mixing effectively depends on the spin and decreases as a function of 
it. Namely, the result of the conhguration mixing directly rehects the selfconsistent change 
of the mean-held deformation as increasing the cranking frequency, which is in contrast to 
the case of ^®'^Er, where the mixing probabilities spread and the calculated B{E2) values by 
the conhguration mixing do not directly rehect the selfconsistent change of the mean-held 
deformation. 

4 Summary and Discussion 

In the present work we propose a method to reliably calculate the rotational band by 
employing the fully microscopic angular momentum projection technique combined with 
the cranking model. The method is based on the conhguration mixing after the projection. 
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in the sense of the GCM, employing several cranked HFB states with different rotational 
frequencies; the original idea was suggested by Peierls-Thouless [l^. By applying the method 
to a few realistic examples, it has been found that the necessary number of cranked HFB 
states employed in the conhguration mixing is rather small, something like four to hve in 
the examples studied. The result is essentially independent of the chosen set of cranking 
frequencies and the set of equidistant frequencies in a suitably chosen range has been mainly 
utilized. 

One of the examples of application is the ground state rotational band of a typical rare- 
earth nucleus ^®'^Er; with the Gogny DIS interaction, a reasonably good agreement with the 
experimental data has been obtained for both the spectrum and B{E2) in the whole spin 
range, 0 < / < 22. Especially, the increase of the moment of inertia as a function of spin is 
almost reproduced as a result of the conhguration mixing, while the calculated inertias by the 
simple projection from one cranked HFB state with various frequencies are almost constant. 
Note that there are no adjustable parameters in this calculation. Another example is the 
superdeformed rotational band in ^®^Dy; the fairly constant moment of inertia in a wide spin 
range appears as a result of the conhguration mixing in accordance with the experimental 
data, although the agreement is not perfect because the selfconsistent deformation is slightly 
overestimated. The calculated inertias by the simple projection from one cranked HFB state 
with various frequencies decrease as spin increases in this case; the one with the non-cranked 
HFB state is almost constant but its absolute value is only about 70% of the calculated inertia 
by the projected conhguration mixing. The feature that the calculated inertias decrease is 
also observed in the ground state rotational band of the unstable neutron-skin nucleus "^^Mg 
if calculated by the simple projection from one cranked HFB state. In these two cases, the 
superdeformed band of ^®^Dy and the ground state band of "^^Mg, the pairing correlation 
vanishes for both neutrons and protons; the reduction of inertia may be related to this 
feature. 

In all three examples, the semiclassical cranking moment of inertia is compared with the 
calculated inertia by the projected conhguration mixing. In the superdeformed band of ^^^Dy 
and the ground state band of ^®Mg the two inertias agrees quite well. In the ground state 
band of ^®'^Er, however, the cranking inertia more rapidly increases than the one calculated 
by the projected conhguration mixing, and their diherence reaches about 20% at higher spins. 
Thus only the results of the projected multi-cranked conhguration mixing are in accordance 
with the ex per imental data and clearly indicate the importance of conhguration mixing. 


In Ref. 


13| it was discussed that the superposition with respect to both the position 


and velocity as generating coordinates is necessary in order to reproduce the correct inertial 
mass for the center-of-mass motion; the so-called “double projection” [^. The prescription 
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of Eq. 


is its natural extension to the collective rotational motion 


131. Although the 


inertial mass for rotation is not known a priori, the result of the present work shows that the 
inclusion of the effect of cranking, i.e., superposing the wave functions with different angular 
velocities, is very important. In fact, as it was already stressed in Ref. 1^, the moment of 
inertia calculated by the simple projection from the non-cranked (time-even) HFB state is 
too small. The calculated inertia considerably increases if the projection is performed from a 
cranked HFB state with small cranking frequency; it is the crudest approximation to Eq. (|3|), 
i.e., the frequency integration is approximated by only one point. When the number of points 
is increased the result of the projected conhguration mixing converges independently of the 
chosen set of cranking frequencies, which is natural because the procedure is a discrete 
approximation of the frequency integration in Eq. (E]). What is surprising, however, is that 
the necessary number of cranked HFB states is rather small and the procedure is tractable, 
which is an important consequence of the present work. 

It may be worth mentioning that the behavior of convergence when increasing the number 
of randomly chosen conhgurations for ^®‘^Er is considerably different from those for '^*^Mg and 
i52Dy. Moreover, the necessary number of conhgurations for convergence is largest in the case 
of ^®'^Er. Although we do not understand the dehnite reason for this behavior, it may be 
most probably related to the presence of pairing correlations. 

We have shown that the superposition with respect to the cranking frequency is a key 
for a reliable description of the collective rotation by employing the angular momentum 


projection method. The Projected Shell Model [3| is also based on the projection and the 
conhguration mixing, but does not rely on the cranking model: It utilizes much larger but 
simple conhgurations in a sense of the shell model, i.e., two, four, and, if necessary, higher 
numbers of quasiparticle excited conhgurations associated with the deformed and supercon¬ 
ducting mean-held. We believe that the cranking model is an efficient way to incorporate 
the multi-quasiparticle excitations suitable for describing the collective rotation. To under¬ 
stand the mutual relationship between the Projected Shell Model and the multi-cranked 
conhguration mixing in the present work may be an interesting future problem. 
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